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Unification of the physics of nucleons and nuclei 
O ■ Ulf-G. MeiBner*^ 

o ■ 

(N ■ ^Universitat Bonn, HISKP (Th), Nufiallee 12-14, D-53115 Bonn, Germany 

I outline an ambitious program which aims to achieve a unified description of nucleon and 

■ nuclear properties based on one chiral effective field theory. 
(N ■ 

1. WHAT AND WHY ? 

Tremendous progress has been made in calculating the properties of few-nucleon systems 
^ ' to a very high level of accuracy, as witnessed by many talks during this conference. 

■ In particular, precise calculations of few-body systems are needed since light nuclei are 
O . (almost) the only laboratory for extracting the elusive properties of the neutron. However, 
Q ! the highly successful so-called standard approach has reached certain limitations, for 

example two- and three-nucleon forces are not calculated consistently with each other 
and the physics underlying the nucleon structure is treated very differently from the 

■ one leading to the nuclear bound and scattering states. This induces some theoretical 
2 ■ uncertainty, which is difficult to quantify. For these and other reasons, we really would 

like to have a description of nucleons and nuclei from one field theory. Such a scheme must 
be rooted in the symmetries of QCD and, being based on field theory, should allow for 
^ I a straightforward implementation of gauge and other symmetries. Furthermore, we wish 

■ to have at our disposal some small parameter (s), leading to a systematic and controlled 
expansion. Consequently, theoretical error bars take over the often (mis)used x^/dof = 
1 criterion. Furthermore, such a scheme should be extendable to cope with inelastic 
processes and the onset of relativity. This can be formulated as a new paradigm: To 
calculate nucleon as well as nuclear properties and dynamics, use quantum field theory and 
organize the theory in a systematic expansion based on (a) small parameter (s). Needless 
to say that the working tool will be an appropriately tailored effective field theory (EFT). 
Now one might ask the question whether such a program is not too ambitious from the 
beginning? I do not believe so, chiefiy because QCD itself offers a rich and intricate pattern 
of symmetries and their violations to achieve this goal. More precisely, the spontaneous, 
explicit and anomalous (chiral) symmetry breaking that QCD is supposed to undergo 
can be analyzed in terms of EFT. Furthermore, the recent developments in formulating 
a nuclear EFT allow for the required link between nucleon and nuclear physics, as I 
will show in the following sections. In this short write-up, I can only give a flavor of 
these developments but I hope that the reader will become sufficiently curious. For 
more introductory material and reviews on these various intertwined topics I refer to the 
"handbook of QCD" [1]. 
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2. HOW ? 



For the hght quark flavors, the QCD Hamiltonian possesses an approximate chiral symme- 
try. In the chiral hmit of vanishing quark masses, this symmetry is exact. It is, however, 
not reahzed in the ground state or the particle spectrum. Rather, the chiral symmetry is 
broken down to its vectorial subgroup with the appearance of three pseudoscalar Gold- 
stone bosons, the pions (I restrict myself here to the two flavor case). If one considers 
reactions with at most one nucleon in the asymptotic in- and the out-states, a perturba- 
tive expansion in terms of small momenta and quark masses is possible. Denoting by Q 
a generic small parameter, any S-matrix element or transition current can be written as 



where A is a renormalization scale, Qi denotes coupling {low-energy) constants of the 
contact interactions appearing at each order, and the function / is of order one. The im- 
portant observation is that the index v is bounded from below because of chiral symmetry. 
Thus, provided Q is parametrically small, this defines a perturbative expansion in terms 
of tree and (pion) loop graphs (see also the discussion on the structure of the nucleon 
below). For systems with baryon number larger than one, an additional non-perturbative 
rcsummation is necessary to generate large scattering lengths and shallow nuclear bound 
states, as described in the talks by Beane [2] and Epelbaum [3] at this conference. Next, 
I briefly discuss the structure of the effective chiral Lagrangian, which decomposes into 
the purely mesonic, the meson-nucleon and the nucleonic sectors. It has the generic form 

+ ..., (2) 

where the ellipsis stands for terms with six or more nuclcons fields. Each of these terms 
has a chiral expansion, for example Cj^n = + C^^lf + C^^j^ + . . . , where the superscript 
denotes the chiral dimension. At each order, new and more contact interactions accompa- 
nied by so-called low-energy constants (LECs) appear. While the number of such terms 
quickly rises as one goes to higher orders, the number of operators contributing to a cer- 
tain process stays manageable. Consider e.g. pion-nucleon scattering, analyzed in terms 
of the chiral effective pion-nucleon Lagrangian. We have to determine 4, 5, and 4 LECs at 
dimension two, three and four, respectively, from the total of 7, 23 and 118 terms at these 
orders [4]. Furthermore, the numerical values for these LECs are understood in many 
cases, rooted in the so-called resonance saturation hypothesis (that is integrating out 
heavier degrees of freedom and estimating LECs from the couplings and masses of these 
heavy states), first investigated in the meson sector [5,6]. As an example, take the finite 

(2) 

dimension two coefficient C4 from C)^l^. From the analysis of various vrA^ scattering PWAs, 
one finds C4 = 3.36. . .3.73GeV~^ In terms of resonance saturation [7] (see Fig. 1), we 
get c4 = cf + cf + 4 = (1.9 + 0.1 + 1.6) GeV-^ = 3.6 GeV-^ where A, R and p denote 
the Delta, the Roper resonance and the p meson, in order. We thus see that not only the 
numerical value of this LEG can be understood but also it comes out somewhat larger 
than the value based on naive dimensional analysis (NDA), |c4°^| ~ 1/m ~ 1 GeV^^ 
(with m the nucleon mass). For a similar investigation of the LECs in Cnn (four-nucleon 
operators) see [8]. 





(1) 
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Figure 1. Resonance saturation for a di- 
mension two operator (filled circle). Solid 

and dashed lines stand for nucleons and pi- 
ous, A'^* (M) denotes a baryon (meson) res- 
onance. 

3. FIRST EXAMPLE: NUCLEON AND DEUTERON ELECTROMAGNETIC 
FORM FACTORS 

As a first example how such a unified approach works, we consider elastic electron scat- 
tering off the nucleon and off deuterium. First, we discuss the nucleon case. The nucleon 
matrix element of the conserved vector current is parameterized in terms of two form 
factors Fi^2 (or, equivalently, the Sachs form factors, Ge,m)- These form factors are also 
a good example of how the complicated structure of the nucleon, which is certainly an 
extended object, arises from the structureless spin-1/2 fields in the chiral Lagrangian. 
There arc two distinct contributions. The first finite size effect is related to the already 
mentioned fact that the effective Lagrangian contains all terms allowed by PCT transfor- 
mations and chiral symmetry. This leads to a string of terms of higher order couplings in 
the various sectors. These couplings parameterize all short distance physics contributing 
to the internal structure of the baryons. 




(1) (2) 



(3) 



(4) 



Figure 2. Chiral expansion of a nucleon electromagnetic form factor. The lowest order 
graph (1) with dimension one insertions only gives the charge of the baryon. Tree graphs 
with insertions from the dimension two (2), three (3) and four (4) effective Lagrangian 
are depicted by the solid circle, square and sextanglc, in order. Pion loop graphs start at 
third order (3), and fourth order loop graphs (4) have exactly one dimension two insertion. 
Only one typical loop graph at the orders considered is shown. Solid, dashed and wiggly 
lines denote nucleons, pions and photons, respectively. 



The other finite size effect of the nucleon only starts to get generated at the one-loop level, 
as in every quantum field theory short-lived fluctuations can occur. In baryon CHPT a 
nucleon typically emits a pion, this energetically forbidden ttN intermediate state lives 
for a short while and then the pion is reabsorbed by the nucleon, in accordance with 
the uncertainty principle. This mechanism is responsible for the venerable old idea of 
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Figure 3. Left panel: Neutron electric form factor G'%{Q'^) at third and fourth order in 
comparison to the data, see [10]. Right panel: Deuteron structure function B versus the 
photon virtuality Q at LO, NLO and NNLO in comparison to the data, see [12]. 

the "pion cloud" of the nucleon, which in CHPT can be put on the firm ground of field 
theoretical principles. This is graphically depicted and further explained in Fig. 2. It is, 
however, important to stress that the relative strength of these contributions in general 
depends on the choice of renormalization scale, for a detailed discussion I refer to [9]. In 
the left panel of Fig. 3 I show the result for the neutron electric form factor based on a 
fourth order calculation employing infrared regularization [10] in comparison to modern 
data extracted from polarization experiments at NIKHEF, MAMI and MIT-Bates. There 
is one LEG, which can be fixed from the neutron charge radius measured in thermal 
neutron-atom scattering as indicated by the dashed hne in the figure. As discussed in 
detail in [10], to achieve an equally accurate description of the other three form factors 
for Q"^ < 0.4 GeV^, vector mesons must be included. For these 3-point functions, that 
can be done systematically and with no new tuneable parameters. Space forbids to show 
the resulting form factors, the interested reader should consult Ref. [10]. Let us now turn 
to electron scattering off the deuteron. The central object is the (unpolarized) scattering 
cross section, which is given in terms of two structure functions. 



with = —Q^ < the invariant momentum transfer squared and 6 is the scattering 
angle in the centre-of-mass frame. Furthermore, we have separated the QED (Mott) cross 
section. These structure functions are subject to a chiral expansion, and consist of essen- 
tially two types of contributions. The first type of terms comprise the so-called impulse 
approximation contribution, which is nothing but the embedding of the single nucleon 
form factors within the deuteron, as indicated by the dashed boxes in the upper row of 
Fig. 4. In addition, there are the so-called three-body corrections (meson exchange cur- 
rents), which are parameter-free up to NLO and have one four-nucleon-photon operator 
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at NNLO [11,12]. The corresponding LEG can be fixed at zero momentum transfer from 
the deuteron magnetic moment. In the right panel of Fig. 3 the structure function B{Q) 
is shown at LO, NLO and NNLO. Visibly, the description of the data improves at higher 
orders and, furthermore, the corrections become smaller when going from LO to NLO 
to NNLO. As stressed in [12], the limiting factor in the accuracy to be achieved for the 
deuteron structure functions stems from the description of the single nucleon form factors. 
For a somewhat different approach, see [13]. 



LO 



NLO 



NNLO 



Figure 4. Chiral expansion of ed scat- 
tering. The upper row refers to the 
single scattering contributions (im- 
pulse approximation), the lower row 
to the three-body corrections (me- 
son exchange currents). The triangle 
symbolizes the deuteron wave func- 
tion, for further notation see Fig. 2. 
Only one representative diagram is 
shown at each order. 

While not very spectacular, I think that this example nicely shows how indeed the nucleon 
and nuclear properties can arise from one effective field theory and it will be of utmost 
importance to extend this program to the investigation of the electromagnetic structure 
of light nuclei. 

4. SECOND EXAMPLE: PION-NUCLEON AND PION-DEUTERON 
SCATTERING 



Another interesting example concerns pion-nucleon scattering at zero energy, i.e. the 
S-wave scattering lengths a+ (isoscalar) and a~ (isovector). Already Weinberg pointed 
out many decades ago that these observables allow for a crucial test of chiral symmetry 
breaking. In the mean time, the chiral expansion of these quantities has been worked out 
to fourth order [14,15] (for brevity, we do not display the fourth order contribution to a"^ 
but refer the reader to [15]) 
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(4) 
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where ^ = M^/m, ~ 1/7 is the small threshold parameter, A = —Aci + 2c2 + 2c3 and 
D = di + d2 + + 2(^5 are combinations of dimension two and three LECs, respectively. 
These can be determined by a fit to the existing PWAs at finite energy, so that one 
can predict the scattering lengths and the phases at higher energies. This leads to the 
results collected in Table 1 [15]. Wc sec that for a~ the convergence is good and the 
value has a small theoretical uncertainty (the low value from Fit 2 should be excluded 
as discussed in detail in [16]). On the other hand, the convergence for a"*" is moderate 
and the theoretical error quite large. These results are consistent with the precise pionic 
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Table 1 

Convergence of the S-wave scattering lengths for CHPT fits to the Karlsruhe (Fit 1), 
Ziirich (Fit 2), and VPI (Fit 3) phase-shift analyses of vrA^ data. 0{p") means that all 
terms up-to-and-including order n were included. Units are M~^. 





0{p) 


0{p') 


0{p') 


0{p') 


Fit 1 
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-0.0096 
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0.0790 
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0.0903 


a- Fit 2 


0.0790 


0.0790 


0.0772 


0.0771 


Fit 3 


0.0790 


0.0790 


0.0870 


0.0867 



hydrogen data obtained at PSI, which again do not allow to pin down a'^ very precisely 
(see the discussion in [17]). It is also important to stress that only with the recent work of 
the Bern group [18] we really have a clear-cut analysis of the isospin breaking effects that 
are so important to properly deduce the strong contribution to the energy shift of the 7r~p 
atom (note also that the conventional methods of subtracting the electromagnetic effects 
from low energy nN scattering data has been put in doubt recently [19]). On the other 
hand, the deuteron is an isoscalar target and thus to leading order entirely sensitive to a'^. 
Beyond leading order, rescattering and other type of three-body corrections play a role, 
and these are also sensitive to the isovector scattering length. This problem has already 
been addressed to third order in Refs. [20,21] and to fourth order in [17]. Surprisingly, 
at fourth order there is only one novel type of correction related to boosting the single 
scattering amplitudes (cf. the contribution Oboost in Fig- 5) , all three-body corrections at 
this order vanish (for details see [17]). However, as it is known since long, the rescattering 
contribution (the graph ag^j^^^jy in Fig. 5) is much larger than the other terms at this 
order (the graphs ai'^body i^ Fig. 5), see e.g. [22]. On the other hand, in the standard 
CHPT counting, numerically large contributions like e.g. multiple rescattering (the graph 
^3-body i^ Fig. 5) are suppressed. Consequently, one has to deal in a systematic fashion 
with this issue. This was done in [23], were it was shown that the rescattering type graphs 
are IR enhanced due to the small binding momentum in the deuteron, 

7 = y/ruD Ed/2 = 45 MeV < ~ 140 MeV , (6) 

with mo {Ed) the deuteron mass (binding energy). More precisely, whenever the ex- 
changed pions can go on (close to) their mass shell, as it is the case for zero energy 
scattering, one has an additional scale to cope with. This can be accounted for by a 
modified power counting, where one counts typical momenta in the deuteron as order p^, 
with p the generic small parameter of CHPT. Amphtudes with powers of are then 
suppressed and, denoting by Q the small parameter in the modified counting, one obtains 
an ordering of the various contributions that is consistent with the numerical evaluation 
of the different mechanisms, as depicted in Fig. 5. Within this new scheme, there are 
only a few more diagrams to calculate. At 0{Q^), one has graphs hke 2a with either two 
dimension two or one dimension three (contact term and loop) insertion. At 0{Q^), the 
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Figure 5. Various contributions to the real part of the 7rd scattering length as explained 
in the text. The filled circle denotes an insertion from C^^lf- 



only new contribution stems from a|_body) ^^S- ^- With that, the final formula for the 
real part of 7rd scattering length, valid to 0{Q'^) in the modified power counting, is 



Re a 



wf 



+ + {a+f - 2{a-)\a+ - a-) 



1 / 1 



- 47r \ |g 



wf 



+ aboost + 0{Q^) . (7) 



where we have subsumed higher order effects into the nN scattering lengths and (. . .)wf 
denotes a momentum space deuteron matrix element. Notice that since local two-pion- 
four-nucleon operators are not enhanced in the modified power counting they appear at 
the same order as in baryon CHPT, namely fifth order, and do not affect the fourth-order 
result given in Eq. (7). The 7r~d atomic level shift has been measured to high precision 
at PSI [24], yielding the complex-valued scattering length a-„d = [— 0.0261(±0.0005) + 
i 0.0063(±0.0007) ]M~^. One thus can analyze Eq. (7) using the real part of this number 
in the (a+ — a~)-plane, as shown by the red bands in Fig. 6. This theoretical uncertainty is 
due to the NLO EFT wave functions employed. If one uses the most recent and systematic 
determination of the pionic hydrogen level shift [18] (the green band in Fig. 6) together 
with the empirical constraint from the hydrogen width shown by the grey band in Fig. 6 
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(which will be affected by further experimental [25] as well as theoretical corrections), one 
can determine the S-wave vrA^ scattering lengths 

a- = [0.0936(±0.0011) JM"^ ; a+ = [-0.0029(±0.0009) ]M^^ . (8) 

If one instead were to use the blue band for the hydrogen shift (based on incomplete 
potential model calculations), these numbers would read a~ = [0.0917(±0.0013) ]M~^ 
and a+ = [-0.0034(±0.0008) ]M-\ The theoretical error bars from both determinations 
should be considered somewhat optimistic given the neglect of isospin violation in [17]. 
For a comparison to other (more model-dependent) determinations and a discussion of 
the imaginary part of the nd scattering length the reader is referred to [17]. For a first 
attempt towards a systematic calculation of pion-^He scattering see [26]. 
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Figure 6. Determination of the S- 
wave ttN scattering lengths from pi- 
cnic hydrogen and deuterium as ex- 
plained in the text. The heavy dot 
refers to the current algebra predic- 
tion and the black box to the CHPT 
determination from scattering data 
[16,15]. 



5. SUMMARY AND OUTLOOK 

I have outlined a program that allows to calculate nucleon and nuclear properties from 
one chiral effective field theory, based on a systematic power counting, embodying the 
symmetries of QCD and allowing for the estimation of theoretical uncertainties. As ex- 
amples, I have discussed electron scattering off nucleons and deuterium and also pion- 
nucleon and pion-deuteron scattering lengths, both cases displaying an intricate interplay 
between chiral nucleon and chiral nuclear dynamics. Other examples, which have gained 
some popularity over the last decade, are pion photo/electroproduction (see e.g. the dis- 
cussion in [27] as weU as [28] and [29]) and Compton scattering off nucleons [30] and light 
nuclei [31], to name just two. It remains to be seen how far this program can be carried 
out in terms of mass number A (the extension to systems with A = 3,4 is under way), 
photon virtuality or pion energy u. One also might ask at what point other active de- 
grees of freedom (the delta, vector mesons, . . .) have to be taken into account explicitely? 
Many pioneering calculations have been done but much more work is required to really 
develop this scheme into a true working tool for precision nuclear, hadronic and few-body 
physics. 
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